Astronomy & Astrophysics manuscript no. sn2004dj 'dust'fin4'lang 


©ESO 2010 


December 23, 2010 





Dust formation in the ejecta of the Type II-P supernova 2004dj 

T. Szalai 1 , J. Vinkd 1 , Z. Balog 2 , A. Gaspar 3 , M. Block 3 , and L L. Kiss 4 - 5 

Department of Optics and Quantum Electronics, University of Szeged, Dom ter 9., Szeged H-6720, Hungary 
e-mail: szaszi@titan . physx . u-szeged . hu ; vinko@titan . physx . u-szeged . hu 
Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-691 17 Heidelberg, Germany 
Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 

Konkoly Observatory of the Hungarian Academy of Sciences, H-1525 Budapest, P.O. Box 67, Hungary 
Sydney Institute for Astronomy, School of Physics A28, University of Sydney, NSW 2006, Australia 



o 

(N 



Received ..; accepted .. 

ABSTRACT 



(N 



O . Aims. Core-collapse supernovae (CC SNe), especially Type II-Plateau ones, are thought to be important contributors to cosmic dust 

production. SN 2004dj, one of the closest and brightest SN since 1987A, offered a good opportunity to examine dust-formation 
processes. To find signs of newly formed dust, we analyze all available mid-infrared (MIR) archival data from the Spitzer Space 
Telescope. 

Methods. We re-reduced and analyzed data from IRAC, MIPS, and IRS instruments obtained between +98 and +1381 days after 
explosion and generated light curves and spectra for each epoch. Observed spectral energy distributions are fitted with both analytic 
and numerical models, using the radiative-transfer code MOCASSIN for the latter ones. We also use imaging polarimetric data 
f^j , obtained at +425 days by the Hubble Space Telescope. 

Results. We present convincing evidence of dust formation in the ejecta of SN 2004dj from MIR light curves and spectra. Significant 
MIR excess flux is detected in all bands between 3.6 and 24 pm. In the optical, a ~0.8 % polarization is also detected at a 2-sigma 
level, which exceeds the interstellar polarization in that direction. Our analysis shows that the freshly-formed dust around SN 2004dj 
can be modeled assuming a nearly spherical shell that contains amorphous carbon grains, which cool from ~700 K to ~400 K between 
+267 and +1246 days. Persistent excess flux is found above 10 pm, which is explained by a cold ( — 115 K) dust component. If this cold 
dust is of circumstellar origin, it is likely to be condensed in a cool, dense shell between the forward and reverse shocks. Pre-existing 
circumstellar dust is less likely, but cannot be ruled out. An upper limit of ~8 x 10~ 4 M G is derived for the dust mass, which is similar 
C/3 , to previously published values for other dust-producing SNe. 

a , 
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^ 1 1. Introduction in distant galaxies. Valiante et al. (2009) showed that differ- 

' ent star-formation histories of galaxies should be taken into ac- 

O ! Core-collapse supernovae (CC SNe) are usually interpreted as count to avoid underestimating the contribution of AGB stars to 

. the endpoints in the life cycle of massive stars (M > 8 M ; e.g. the whole dust production. Another possibility is the condensa- 

• ; Woosley et al. 2002). While it is widely accepted that these high- t ion of dust grains in quasar winds (Elvis et al. 2002), which 

, energy events have significant effects on the birth and evolution has also been supported by observations (Markwick-Kemper et 

q ' of other stars, there has been a long-standing argument whether a i. 2007). On the other hand, there are many high-z galaxies 

; they also play an important role in the formation of interstellar wne re SN-explosions are the only viable explanation for the 

!* ■ dust. observed amount of dust (Maiolino et al. 2004; Stratta et al. 

Theories on the dust-production of CC SNe have been 2007; Michalowski et al. 2010). Detailed studies by Matsuura et 

^ . around for over 40 years (Cernuschi et al. 1967; Hoyle & al. (2009) suggest a 'missing dust-mass problem' in the Large 

5_i ' Wickramasinghe 1970). These early hypotheses were later sup- Magellanic Cloud (LMC), similar to the situation in distant 

ported by studies of isotopic anomalies in meteorites (Clayton galaxies described above. 

1979; Clayton etal. 1997; Clayton &Nittler 2004). Additionally, The two main problems with the observable dust content 

several papers were published about the unexpectedly large dust around SNe are the amount and the origin. Models of Kozasa 

content of high-redshift galaxies (Pei et al. 1991; Pettini et al. e t al. (1989), Todini & Ferrara (2001), and Nozawa et al. 

1997; Bertoldi et al. 2003), which suggested that CC SNe could (2003) predict 0.1-1 M of newly formed dust after a CC SN 

be the main sources of interstellar dust in the early (and maybe event. More recently, Bianchi & Schneider (2007), Kozasa et al. 

in the present) Universe (Todini & Ferrara 2001; Nozawa et al. (2009), and Silvia et al. (2010) have found similar values via 

2003; Morgan & Edmunds 2003). The average lifetimes of CC numerical modeling of the survival of dust grains. They have 

SNe are short enough to produce dust at high redshifts, i.e. dur- found that the dust amount depends on the type of SN and also 

ing ~1 Gyr, thus these SNe are better candidates for explaining the density of the local ISM, which must be taken into account 

the presence of early dust than other possible sources, like AGB when estimating the contribution of CC SNe to the dust content 

stars (e.g. Dwek et al. 2007). in high-redshift galaxies. 

Although these theories have not been disclaimed so far, Despite the nearly concordant results of different models, di- 

there are several other possible mechanisms of dust formation rect observations have not confirmed the massive dust produc- 
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tion by CC SNe yet (although a detailed analysis was possible 
only in a few cases). The first evidence for dust condensation was 
observed in SN 1987A (e.g. Danziger et al. 1989, 1991; Lucy 
et al. 1989, 1991; Roche et al. 1993 and Wooden et al. 1993; 
revisited by Ercolano et al. 2007). The observed signs of dust 
formation were i) a strong decrease of optical fluxes around 500 
days after explosion, if) an increase of mid-infrared (MIR) fluxes 
at the same time, and Hi) increasing blueshift and asymmetry of 
the optical emission lines because of an attenuation of the back 
side of the ejecta by the freshly synthesized dust. The mass of 
newly condensed dust was estimated to be ~ 10 -4 M Q . Similar 
effects were also detected in the ejecta of Type II-P SN 1999em 
and Type Ib/c SN 19901 (Elmhamdi et al. 2003, 2004). 

The launch of the Spitzer Space Telescope (hereafter Spitzer) 
in 2003 and AKARI in 2006 provided additional opportunities to 
observe dust around CC SNe by following the evolution of MIR 
light curves and spectra. Evidence for newly condensed dust 
were obtained for SN 2003gd (Sugerman et al. 2006; Meikle 
et al. 2007), 2004et (Kotak et al. 2009), 2007od (Andrews et al. 
2010), and Type Ib/c SN 2006jc (Nozawa et al. 2008; Mattila et 
al. 2008; Tominaga et al. 2008 and Sakon et al. 2009). In every 
case, the estimated mass of recently formed dust was between 
10~ 5 - 10~ 3 M G . Note that for SN 2003gd Sugerman et al. (2006) 
calculated a value of 0.02 M , but their result was questioned by 
Meikle et al. (2007). These values are significantly lower than 
the theoretically predicted ones and tend not to support an SN 
origin for observable amounts of dust in the local and distant 
universe. 

There are other ways to discover dust around SNe. As it 
was found in some CC SNe (e.g. SN 1998S, Pozzo et al. 2004; 
SN 2005ip, Smith et al. 2009 and Fox et al. 2009; SN 2006jc, 
Smith et al. 2008a; SN 2006tf, Smith et al. 2008b; SN 2007od, 
Andrews et al. 2010), dust grains may condense in a cool dense 
shell (CDS) that is generated between the forward and reverse 
shock waves during the interaction of the SN ejecta and the pre- 
existing circumstellar medium (CSM). The CDS may affect both 
the light curves and the spectral line profiles. Another possi- 
bility is the thermal radiation of pre-existing dust in the CSM 
that is re-heated by the SN as an IR-echo (Bode & Evans 1980; 
Dwek 1983, 1985; Sugerman 2003), which can be observed as 
an infrared excess (e.g. SN 1998S, Gerardy et al. 2002, Pozzo 
et al. 2004; SN 2002hh, Barlow et al. 2005, Meikle et al. 2006; 
SN 2006jc, Mattila et al. 2008; SN 2004et, Kotak et al. 2009; 
SN 2008S, Botticella et al. 2009). These results seem to support 
the hypothesis that pre-explosion mass-loss processes of progen- 
itor stars could play a more important role in dust formation than 
the explosions of CC SNe (see also Prieto et al. 2008 and Wesson 
etal. 2010). 

Recent studies of SN remnants (SNRs) could not solve the 
question of the amount of dust produced by SNe. By obtain- 
ing far-IR and sub-mm observations, many groups estimated the 
amount of newly condensed dust grains in Cas A (Dunne et al. 
2003; Krause et al. 2004; Rho et al. 2008). Their results varied 
between 0.02 and 2 M , while the values for Kepler SNR showed 
an even larger discrepancy (0. 1 -3 M e by Morgan et al. 2003, and 
5 x 10~ 4 M by Blair et al. 2007). Stanimirovic et al. (2005) and 
Sandstrom et al. (2009) studied the MIR data of SNR 1E0102.2- 
7219 in the Small Magellanic Cloud (SMC) and found evidence 
for ~ 1 — 3 x 10~ 3 M Q of dust formed in the ejecta (the authors of 
the latter paper also calculated that the amount of cold dust that 
is observable only at longer wavelengths could be even higher 
by two orders of magnitude). 

Estimating the mass of dust around SNe is complicated, and 
the results are strongly model-dependent. Several ideas have 



been proposed to explain the discrepancies between observations 
and theory. Using clumpy grain-density distributions (Sugerman 
et al. 2006, Ercolano et al. 2007) leads to a dust mass in a SN 
ejecta at least one order of magnitude higher than by assuming a 
smooth density distribution. Regarding dust enrichment in high- 
z galaxies, the low dust production rates of observed CC SNe 
might be compensated by using a top-heavy initial mass function 
(IMF), resulting in more SNe per unit stellar mass (Michalowski 
et al. 2010 and references therein). It is also a possibility that SN 
explosions provide only the dust seeds, while the bulk of the dust 
mass is accumulated during grain growth in the ISM (Draine 
2003, 2009; Michalowski et al. 2010 and references therein). It 
has also been noted that estimates of the dust content at high 
redshifts should be considered with caution because of the un- 
certainties in interpreting low signal-to-noise observations (see 
e.g. Zafar etal. 2010). 

The aim of this paper is to study the dust formation in the 
nearby bright Type II-P SN 2004dj. This SN occured in the com- 
pact cluster Sandage-96 (S96) in NGC 2403, and owing to its 
proximity (~ 3.5 Mpc, Vinkd et al. 2006) it has been intensively 
studied since the discovery by K. Itagaki (Nakano et al. 2004, 
Patat et al. 2004). Early observations have been summarized by 
Vinko et al. (2006, hereafter Paper I), while the physical proper- 
ties of S96 have been recently discussed by Vinko et al. (2009, 
hereafter Paper II). The progenitor was identified as a massive 
(12 M Q < M prog < 20 M Q ) star within S96 (Mafz-Apellaniz et 
al. 2004; Wang et al. 2005; Paper II). 

The evolution of SN 2004dj was extensively observed by 
Spitzer. The earliest observations have been presented by Kotak 
et al. (2005). In the following sections we analyze the Spitzer 
and Hubble observations that show signs of dust formation 
(Section [2] and Section [3), then we compare various dust mod- 
els with these observations in Section [4] Finally, we present our 
conclusions in Section|5] 

2. Observations and data reduction 

This section contains the description of the Spitzer- and Hubble 
observations that we use to probe the newly formed dust in 
SN 2004dj. Throughout this paper we adopt JD 2,453,187.0 
(2004-06-30) as the moment of explosion (see Paper I). 

2.1. MIR Photometry with Spitzer 

We collected all public archival Spitzer data on SN 2004dj for 
each channel (3.6, 4.5, 5.8 and 8.0 pm) of the Infrared Array 
Camera (IRAC) and for the 24 pm channel of the Multiband 
Imaging Spectrometer (MIPS). We use all available basic- and 
post-basic calibrated data (BCD and PBCD) obtained from +98 
to +1381 days past explosion (summarized in Table[TJ. 

The IRAC BCD images were processed using the IRACproc 
software (Schuster et al. 2006) to create final mosaics with a 
scale of 0.86"/pixel. Aperture photometry was also carried out 
with this package (which uses built-in IRAB3 scripts) with an 
aperture radius 2" and a background annulus from 2" to 6". 
Aperture corrections of 1.213, 1.234, 1.379, and 1.584 were 
used for channels 1-4, respectively (Table 5.7 in IRAC Data 
Handbook, Reach et al. 2006). We compared our results against 
aperture photometry on PBCD images using the IRAF phot task 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Post-BCD images showing the environment of SN 2004dj obtained with Spitzer on 2004-10-12: IRAC 3.6 micron (top left), 
5.8 micron (top right), 8.0 micron (bottom left), and MIPS 24.0 micron (bottom right). The SN is marked in each panel. The target 
is very faint on the MIPS image, which illustrates the necessity of PSF-photometry (see text for details). 



with the same parameters, and found a reasonably good agree- 
ment. 

Registration, mosaicing and resampling of the MIPS images 
were made by using the DAT software (Engelbracht et al. 2007). 
The source was not detected on the 70 pm images, therefore 
we analyzed only the 24 fim frames. The photometry of the 24 
/im data was complicated, because the source was very faint, 
peaking only slightly above background (see FigJTJ. First, aper- 
ture photometry was performed on the 1.245 "/pixel scale mo- 
saics using an aperture radius of 3.5" and a background annu- 
lus from 6" to 8". An aperture correction of 2.78 was then ap- 
plied (Engelbracht et al. 2007). The choice of such a small aper- 
ture and annulus was motivated by the close, bright H II region, 
which makes the background high and variable on one side of 
the source. 



Second, PSF-photometry was carried out on the 24 yum 
MIPS frames using a custom-made empirical PSF (assem- 
bled from multiple epoch observations of HD173398) using 
IRAF/DAOPHOT. We set the PSF radius as 90 pixels, for 
which no computed aperture correction is available. Therefore, 
the photometric calibration was done by comparing the results 
of the properly calibrated aperture photometry and our PSF- 
photometry for two isolated point-like sources within the SN- 
field that do not suffer from the contamination of a nearby 
bright object. We compared the the sum of the count rates 
on the DAT frames (in DN/s) from both aperture- and PSF- 
photometry, aperture-corrected the numbers from aperture pho- 
tometry and found the following linear relationship betweeen the 
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aperture-corrected count rates (Japan) an d those from our PSF- 
photometry (Jpsf)'- 



fapcor = 1.183(±0.012)x/ PSF + 16.53(±16.91). 



(1) 



The fluxes from PSF-photometry corrected this way were 
then converted to mJy as prescribed for MIPS and found to be 
within the errors of those from aperture photometry. 

We also performed PSF-photometry with the IDP3 software 
(Stobie & Ferro 2006), which agreed well with the results of 
DAOPHOT. As a triple-check, we also performed aperture pho- 
tometry on the PBCD-images and found reasonable agreement 
between the absolute flux levels computed from PSF- and aper- 
ture photometry. Thus, the final 24 pm fluxes were calculated as 
the average of the results from the different methods. The final 
photometry of SN 2004dj is collected in TableQ] The errors rep- 
resent the standard deviation of the data from the three different 
photometry methods. 

Kotak et al. (2005) presented some MIR photometric points 
observed at ~ 1 00 days after explosion. Their reported flux values 
are generally consistent with ours, but there are some differences 
(mainly in channels 4.5 and 24 pm). The reason for these minor 
differences may be that they applied larger aperture radii (3.6" 
and 5.6" for 4.5 and 24 /mi data, respectively) and performed 
only aperture photometry on the PBCD images. 

2.2. MIR Spectroscopy with IRS 

SN 2004dj was observed with the Infrared Spectrograph (IRS) 
onboard Spitzer on seven epochs between October 2004 and 
November 2006, from + 1 15 to +868 days after explosion. These 
observations were made in the IRSStare mode using the Short- 
Low (SL) setup. Table [2] contains the list of spectroscopic data 
downloaded from the Spitzer archive. In addition, several imag- 
ing observations were made with the blue array (16 pm central 
wavelength) of the peak-up imaging (PUI) mode of IRS, which 
are listed in Table Q] 

We analyzed the PBCD-frames containing the spectro- 
scopic data using the SPitzer IRS Custom Extraction software 
(SPIClQ). Sky subtraction and bad pixel removal were per- 
formed using two exposures containing the spectrum at different 
locations, and subtracting them from each other. Order extrac- 
tion, wavelength- and flux-calibration were performed applying 
built-in templates within SPICE. Finally, the spectra from the 
1st, 2nd and 3rd orders were combined into a single spectrum 
with the overlapping edges averaged. In a few cases the sky near 
the order edges was oversubtracted because of excess flux close 
to the source, which resulted in spurious negative flux values in 
the extracted spectrum. These were filtered out using the fluxes 
in the same wavelength region that were extracted from the ad- 
jacent orders. The results were also checked by comparing the 
extracted fluxes with photometry (Table [TJ. Reasonable agree- 
ment between the spectral and photometric fluxes was found for 
all spectra. 

The final calibrated and combined spectra cover the 5.15 - 
14.23 pm wavelength regime with resolving power R ~ 100. 
These are plotted in Figf2] where a small vertical shift is applied 
between the spectra for better visibility. The analysis of the spec- 
tral features and evolution will be presented in Section|3] 

The fluxes of SN 2004dj in the PUI frames were measured 
via aperture photometry by the MOsaicker and Point source 
Extractor software. The results are given in TableQ] 

2 http://ssc.spitzer.caltech.edu/dataanalysistools/tools/spice/ 

3 http://ssc.spitzer.caltech.edu/dataanalysistools/tools/mopex/ 
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Fig. 2. IRS spectra of SN 2004dj in the nebular phase. Line iden- 
tification is based on Kotak et al. (2005, 2006). The evolution of 
the features is discussed in Sect. 13.21 



These broad-band integrated fluxes were only used for checking 
the shape of the SED between the IRAC and MIPS bands, but 
were not included in the detailed analysis in Section [3] 



2.3. Imaging polarimetry with HST 

SN 2004dj and its surrounding cluster, Sandage-96, were imaged 
by HST/ ACS High-Resolution Camera (HRC) on 2005 Aug 28 
(proposal GO-10607, PI: B.E.K. Sugerman), +425 days after ex- 
plosion. Among others, three sets of four drizzled frames were 
recorded through the F435W filter at 0, 60, and 120 degree po- 
sition of the POLUV polarization filter (see Paper II for the 
description and analysis of all observations). Here we consider 
only the polarization measurement of the SN. 

We measured the flux from SN 2004dj with PSF-photometry 
using the DOLPHOT software (Dolphin 2000) the same way as 
described in Paper II. From the fluxes measured at three different 
polarizer angles (Iq, I(,q and /120), the Stokes-parameters (/, Q, 
U) were calculated applying the formulae given by e.g. Sparks 
& Axon (1999). The degree-of -polarization, p was then derived 
as 



P = 



VG 2 + u 2 



(2) 
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Table 1. Spitzer photometry for SN 2004dj. t eX pi = 2,453,187 was adopted as the moment of explosion. 



Flux (1(T 20 erg s~ ] crrT 2 A" 1 ) 



U 1 Date 


A /TTT~\ 

JVIJD — 






IRAC 




TO C TJT TT A/fTTiO 




24D0U0U 


(days) 


3.6 /jm 


4.5 /iin 


j.o /Jin 


8.0 yum 


13. U - lo.j yum 24 /jm 


iaa/i 1 a rvia 
ZUU4- 1 U-U / 


JZSJ.D 


no 


24800(jy) 


i jyuu(,j i ) 




ly40(16) 




iaa/i 1 a aq/j 
ZUU4-lU-Uo 


jzBo.y 




24U0U(jy) 


1 /1AAA/T")\ 

14000(3Z) 


jjoO(Z9) 


1oj0(16) 




2004-10-12 


3290.6 


103 


1 o o r\r\f f i \ 

18300(61) 


1 Q/iAA/1A\ 

13oU0(3U) 


CAAA/1 A\ 

j000(Z9) 


1530(15) 




r-\ f\f\ A 1 f\ 1 >~th 

2004-10-12 


3291.4 


104 










45(3) 


2004-10-14 


3293.0 


106 










50(6) 


oah/i 1 a i f h 

2004-10-16 


3295.0 


108 










31(2) 


2004-1 1-01 


3310.6 


123 


A/'OA/OA\ 

9680(39) 


1 ATAA/ t1\ 

10700(32) 


3960(25) 


1 AAA/ 1 1 \ 

1090(11) 




2004-1 1-06 


3316.2 


129 










55(2) 


2005-03-03 


3432.8 


245 










48(2) 


r-\ f\r\ c r\i r\ An 

2005-03-24" 


3454.4 


267 


2740(22) 


A A 1 A/ 1 "7^ 

4910(1 /) 


1 lOA/ 1 OA 

IZoO(lo) 


551(11) 




Z00j-04-01 


3462.5 


275 










44(2) 


ZUUo-lU-ZU 


"2AA/1 1 

3004. Z 


4 / / 


41 1U(ZZ) 


3ZoO(lZ) 


1 A 1 A/1 1 \ 

1910(Z1) 


10j0(1 i) 




Z00j-1 1-ZZ 


3696.9 


510 










138(1) 


iaaa m tjc 


jolo.4 


A "2 1 
03 1 


zysU(20) 


i/i 'jn/ 1 1 \ 
Z430(l 1 ) 


1 AAA/ 1 A\ 
190U(10) 


111 n/ 1 t\ 
1 1 10(12) 




ZUUO-U4-Z3 


1Q AO A 


ooZ 










1 OO/ 1 \ 

182(1) 


1AAA 1 A 10(/ 

ZUOo-lO-Zo 


4030. / 


o49 


18/0(1S) 


lOlU^V; 


14VU(,1 j J 


n i n/ 1 1\ 

yiy(i2) 




ZUUO-lUol 


AfYXQ A 
4U3V.0 


ojZ 


1ojU(j.3) 


1 SSO/^Ch 
1 JoU^jy ) 


IJl U^O J ^ 


Q/i ^/")c^ 




ZUUO- 11-10 


4UJ4. / 


CA7 
50/ 










1/3(4) 


onnA ii i a<? 
ZUUO- 11-10 


/I AC C A 

4UJJ.U 


OAC 

0O0 










1 r 7 r )t r )\ 


1 AAA 11 A1C 

ZUUo- 1 Z-0 1 


Af\Hf\ O 

40/0.0 


oo3 










66(4) 


nam A/1 Al*? 

ZUU /-U4-UZ 


a i m t 
4193.3 


1 AAA 
1UU0 


1U/U(16) 


Q 1 1( 1 C\\ 
y 1 / 1U; 


you^ ) 


A^n/ 1 i \ 
ojU(1 1 ) 




OAAT A/1 Alt? 

ZUU /-U4-UZ 


a i m a 
4193.3 


1 AAA 
1UU0 


lUaU(lo) 


Q99(' 1 1 A 
VZZ^ 1 L) 


O l\J\ 1 J) 


602(10) 




OAA7 A/1 1 Q<? 

ZUU / -U4- 1 J 


*+ZU3.0 


i ni a 

1U10 










DJ(4) 


2007-10-24 / 


4397.8 


1210 










59(4) 


2007-10-24* 


4397.8 


1210 










58(6) 


2007-1 1-04? 


4408.3 


1221 










117(3) 


2007-1 l-19 e 


4423.8 


1236 


779(15) 


616(8) 


619(11) 


437(10) 




2007-1 1-23' 


4427.6 


1240 


771(15) 


613(8) 


619(11) 


444(10) 




2007-11-24* 


4428.5 


1241 


787(14) 


617(7) 


647(7) 


444(10) 




2007-1 l-29 c 


4433.9 


1246 










62(3) 


2007-12-15* 


4449.9 


1262 










113(2) 


2008-04-07* 


4564.3 


1377 


723(15) 


520(6) 


544(11) 


399(11) 




2008-04- 12^ 


4568.7 


1381 


718(16) 


526(7) 


527(13) 


413(11) 





Notes. 

PID 226 Van Dyk et al. (SONS) 
(W PID 159 Kennicutt et al. (SINGS); 24.0 fim values are MlPSScan data 
(c) PID 20256 Meikle et al. (MISC) 

PID 30292 Meikle et al. (MISC) 
(c) PID 30494 Sugerman et al. (BEKS) 
W PID 40010 Meixner et al. (SEEDS) 
<g) PID 40619 Kotak et al. (MISC) 

This resulted in a measured degree-of-polarization p = 
0.0941 + 0.0029, which is, of course, much higher than the true 
p from SN 2004dj, owing to instrumental polarization (/?,), and 
the polarization from interstellar matter (pis m). 

The instrumental polarization of ACS was studied by Biretta 
et al. (2004). Using three sets of in-orbit calibration data they 
measured the instrumental polarization as p t = 0.086 + 0.002 for 
the ACS/HRC+F435W+POLf/y detector+filter combination 
(see their Table 20). They also give a semi-empirical formula 
describing the relative uncertainty of the degree-of-polarization, 
(Tplp, as a function of the signal-to-noise of the flux measure- 
ment. 

Adopting this calibration, the measured true degree-of- 
polarization, corrected for the instrumental zero-point, is p true = 
0.008 1 +0.0036. The average S/N of the three flux measurements 
for SN 2004dj is ~ 200, from which we derived cr p /p = 0.5588, 



thus, cr p = 0.0045, which agrees well with the error estimate of 
p trm given above. 

Further analysis and discussion of the polarization data are 
presented in Section[3] 



3. Analysis of the observations 

3. 1 . The effect of the surrounding cluster Sandage-96 

The bright, compact cluster S96 surrounding SN 2004dj (Paper 
II) made late-time optical measurements difficult, because the 
light from the cluster stars contribute significantly to the total 
measured fluxes. This was very problematic in the optical, but 
fortunately less severe in the MIR, where the cluster is much 
fainter. Although we have no pre-SN MIR photometry for S96, 
we estimated its MIR flux using the best-fitting cluster models 
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Table 2. IRS observations of SN 2004dj 



UT Date 


MJD - 


t - texpl 


ID 


Proposal ID (PID) 




2,450, 000 


(days) 






2004-10-24 


3302.3 


115 


1-12113152 


226 (SONS, Van Dyk) 


2004-11-16 


3325.8 


139 


rl21 14432 


226 (SONS, Van Dyk) 


2005-03-18 


3447.6 


261 


rl21 13408 


226 (SONS, Van Dyk) 


2005-04-17 


3477.9 


291 


rl21 15456 


226 (SONS, Van Dyk) 


2005-11-22 


3696.8 


510 


rl4458880 


20256 (MISC, Meikle) 


2006-04-23 


3848.4 


661 


1-14465280 


20256 (MISC, Meikle) 


2006-11-16 


4055.1 


868 


r 17969 152 


30292 (MISC, Meikle) 



Table 3. Identified lines in the IRS spectra of SN 2004dj. "+" signs represent positive detection, while "-" symbols mean non- 
detection of the given feature at the specified epochs. 



A (p.m) 


Ion 


+ 115d 


+ 139d 


+261d +291d +510d 


+661d 


+868d 


4.55 


CO 


+ 


+ 


+ + - 






6.6342 


[Nill] 


+ 


+ 


+ + + 


+ 


+ 


6.9852 


[Aril] 






+ + + 






7.4578 


H-Pfa 


+ 


+ 


+ - - 






7.5005 


H-HuyS 


+ 


+ 


+ + + 


+ 




7.5066 


[Nil] 






blend with H-Hu/? 






8.7577 


H7-10 


+ 


+ 


+ + - 






10.521 


[Coll] 


+ 


+ 


+ + + 






11.306 


H7-9 






blend with [Ni I] 






11.308 


[Nil] 


+ 




+ + 






12.368 


H-Hua 


+ 


+ 


+ + + 






12.729 


[Nill] 






blend with [Nell] 






12.813 


[Nell] 


+ 


+ 


+ + + 


+ 


+ 





100000 f 


•< 

cm 


10000 r 




1 erg/s/cn 


1000 - 
100 - 






O 


10 - 


Flux 


1 r 



0.1 



0.01 



0.1 




Table 4. Contribution of the Sandage-96 cluster to the MIR 
fluxes of SN 2004dj 



Wavelength (urn) 



Fig. 3. Best-fitting model SEDs for the host cluster Sandage-96 
from Paper II. Filled circles represent observed data for the clus- 
ter, while open symbols denote the averaged MIR-fluxes from 
the models (see Table|4]l. 



from Paper II. These models are simple, coeval stellar popula- 
tions (SSP), but we considered three different model sets using 
different input physics to test the model dependence of the flux 
predictions. 

Figure [3] illustrates the MIR SEDs for the three SSP mod- 
els that produced the best fit to the observed UV+optical+NIR 
data of S96 (Fig[3] filled symbols; see Paper II for details). 



Wavelength 


Model flux 


Flux error 


Oum) 


(erg s" 


_I cm 


- 2 A- 1 ) 




IRAC 3.6 


5.25 x 10" 


is 


6.2 x 10" 


ii) 


IRAC 4.5 


1.89 x 10" 


is 


3.4 x 10" 


19 


IRAC 5.8 


8.40 x 10' 


19 


1.3 x 10" 


19 


IRAC 8.0 


2.59 x 10" 


19 


0.4 x 10" 


19 


MIPS 24.0 


3.37 x 10" 


21 


0.9 x lO- 


-21 



Fortunately, the three models predicted very similar MIR fluxes 
in the Spitzer bands. Their average values (in cgs units) are col- 
lected in Table H] where the errors illustrate the flux differences 
between the models. 

Evidently the cluster contributes significantly to the mea- 
sured SN fluxes only in the IRAC channels, but the cluster flux 
drops below the observational uncertainties longward of 10 /im. 
Thus, we subsequently corrected the IRAC fluxes for this offset 
using the model fluxes from Table [H but this correction was not 
necessary for the PUI and MIPS data. 

3.2. Spectral features and evolution 

The nebular spectra of SN 2004dj presented in Figf2]are typical 
for Type II-P SNe. They show features and evolution very similar 
to the available MIR spectra of a few SNe, including SN 1987 A 
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Table 5. Measured line fluxes and ionization fraction (x) for 
58 Ni and 56 Co 



Phase 
(days) 


[Coll]" 
10.52 yum 


[Nill] 
6.63 /jm 


[Ni I] 6 
7.51 yum 


x( S8 Ni) 


115 


1.6 


8.47 






139 




17.8 






261 


16.1 


28.4 


2.5 


0.87 


291 


16.8 


27.5 


2.0 


0.89 


510 


8.3 


27.9 


1.86 


0.89 


661 




25.4 


1.88 


0.88 


868 


1.1 


16.7 







Notes. 

(fll xl0~ 15 erg s -1 cnr 2 

{b> corrected for Pace, Hu/? and HI 7-11 contribution (see text) 



(Wooden et al. 1993; Roche et al. 1993), SN 2005af (Kotaket al. 
2006) and SN 2004et (Kotak et al. 2009). The first two spectra of 
SN 2004dj (taken at +1 15 d and +139 d) have been already pre- 
sented and discussed by Kotak el al. (2005). We include these 
spectra in the present study for completeness, but the conclu- 
sions for them are generally the same as those of Kotak et al. 
(2005). 

The spectra consist of permitted emission lines of HI 
(Pfund-, Humphrey- and n=7 series), and forbidden lines of 
[Nil], [Nill], [Coll] and [Aril]. Table [3] summarizes the identi- 
fied lines, based on Wooden et al. (1993) and Kotak et al. (2005). 
The presence of [Ne II], forming a blend with [Ni II] at 12.8 /an, 
was suggested by Kotak et al. (2006) for SN 2005af. This fea- 
ture is also present in SN 2004dj for all phases covered by IRS- 
observations, although the two components cannot be resolved. 

The blue edge of the first two spectra is influenced by the 
red wing of the 1-0 vibrational transition of the CO-molecule, 
similar to SN 1987A, as first pointed out by Kotak et al. (2005). 
CO 1-0 remained detectable in the two subsequent spectra up to 
+291 days. No sign of the SiO fundamental band around 8 fim 
was detected in SN 2004dj, unlike in SNe 1987A, 2005af and 
2004et. Because we present evidence that SN 2004dj showed 
significant dust formation after +400 days, the absence of de- 
tectable SiO can be an important constraint for the chemical 
composition of the newly formed dust grains. 

3.3. Masses of Ni and Co 

Masses of freshly synthesized Ni estimated from observations 
may provide important constraints for SN explosion models (e.g. 
Woosley & Weaver 1995; Thielemann et al. 1996; Chieffi & 
Limongi 2004; Nomoto et al. 2006). Following the methods pre- 
sented in Roche et al. (1993) and Wooden et al. (1993), we de- 
rived the ionization fraction of stable 58 Ni from the measured 
ratio of fluxes of the collisionally-excited lines of [Nill] at 6.63 
yum and [Nil] at 7.51 /mi. It was shown by Wooden et al. (1993) 
and Kotak et al. (2005) that the critical density (at which the 
collisional de-excitation rate equals the rate of spontaneous de- 
cay) for the 6.63 /im line is ~ 1.3 x 10 7 ctrT 3 , while the elec- 
tron density ~ 1 year after explosion is ~ 10 8 ctrT 3 , one order 
of magnitude higher, making LTE-conditions valid. Under these 
conditions the total number of excited atoms at the upper level u 
of the given transition can be expressed as 



N u = 47rD 2 F u ,A/hcA u i, 



(3) 



the transition and A u \ is the transition probability (in cm -1 ). Note 
that this expression is valid for optically thin lines, but the [Ni II] 
6.63 yum feature was probably optically thick between 100-660 
days (Wooden et al., 1993). Thus, the derived quantities for this 
line are actually lower limits. The ratio of the total number of 
neutral Ni° and ionized Ni + can be obtained from the ratio of 
measured line fluxes via Boltzmann- and Saha-equations. 

It is important that the 7.51 fim line in the IRS spectra is an 
unresolved blend of [Nil] with Paa, HuyS and HI 7-11, which 
makes the direct measurement of the [Ni I] feature difficult. This 
was already noted by Kotak et al. (2005), but in a subsequent 
paper Kotak et al. (2006) measured the flux of this line directly 
and derived the ionization fraction x = Ni + / (Ni° + Ni + ) ~ 
0.4 for SN 2004dj at +261 days. This is roughly a factor of 2 
lower than what was measured in SN 1987 A at that epoch (x ~ 
0.8, Wooden et al. 1993). In order to clarify this, we estimated 
the line flux of the [Ni I] feature assuming that the H I features 
have the same ratio to the total line flux as was measured by 
Wooden et al. (1993) for SN 1987 A in higher resolution spectra. 
The contribution from H I was found to be as high as 90 percent 
at +260 days, which decreased to 75 percent at +450 days and 
less than 10 percent at +660 days. We used these numbers to 
estimate the [Ni I] fluxes from the measured total line fluxes (the 
[Nill] 6.63 fim feature was not affected by this complication). 
The results are collected in Table [5] Using these line fluxes and 
assuming T e = 3000 K electron temperature (Wooden et al., 
1993), the ionization fraction was computed as 



3' 



/V(Ni + ) 



F(6.63)6.63 g 7 .5iA 7 .5i z + (T e ) 
F(7.51)7.5U 6 . 63 A 6 . 6 3 z°(T e ) 



yv(Ni°) 

exp[(£ 6 .63 - E r5 i)/kT e ], 



(4) 



where D is the distance to the SN (3.5 Mpc, Paper I), F u i is 
the measured line flux (in erg s _1 cm -2 ), A is the wavelength of 



then applying x = y/(l + y), where g is the statistical weight, 
z + '°(T) are the partition functions for ionized and neutral Ni, re- 
spectively, and E is the excitation potential for the upper level of 
the given transition. The values of the atomic parameters were 
adopted from the NIST Atomic Spectral DatabaseQ. The derived 
ionization fractions are shown in the last column of Table [5] 
These values are very close to those obtained by Wooden et al. 
(1993) for SN 1987A, but significantly higher than the results 
by Kotak et al. (2006). It is probable that Kotak et al. overes- 
timated the contribution of [Nil] to the 7.51 yum feature, and 
therefore underestimated the ionization fraction for Ni. Our new 
result suggests that the Ni ionization and probably other physical 
conditions in the ejecta of SN 2004dj were quite similar to those 
in SN 1987A. 

We also derived the total mass of 58 Ni by applying Eq.3 to 
the 6.63 fim feature and assuming LTE. This resulted in M( 58 Ni) 
~ 5 x 10~ 4 M , a factor of 2 higher than Kotak et al. (2005) 
obtained from the first two spectra. Note that the mass of the 
radioactive 56 Ni synthesized during the explosion was about 2 x 
10~ 2 M (Paper I), about two orders of magnitude higher (see 
FigU. 

A similar analysis could not be completed for Co, because 
neither the neutral [Col] forbidden lines between 3.0 and 3.75 
fim (Wooden et al. 1993), nor the [Col] 12.25// feature (Roche 
et al. 1993) were detectable in the low -resolution IRS spec- 
tra. Nevertheless, we derived masses of singly ionized Co via 
the strength of the [Coll] 10.52// transition. This feature was 
weak in the +115 and +139 day spectra, but became one of the 
strongest features in the +261 and +291 day spectra. At +510 



http://www.nist.gov/physlab/data/asd.cfm 
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Fig. 4. Ni and Co masses in M calculated from the strength 
of forbidden emission lines via Eq.3. Open symbols: Ni-masses 
from the [Nill] 6.63// feature; filled symbols: Co-masses from 
the [Coll] 10.52// feature. Lines show the expected amount of 
Co predicted from the radioactive decay of 0.02 M radioactive 
56 Ni synthesized in the explosion. Continuous line: only 56 Co; 
dashed line: 56 Co plus 57 Co with solar abundance ratio; dotted 
line: 56 Co plus 57 Co assuming twice solar abundance ratio. 
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Fig. 5. Temporal evolution of the 10 fim continuum flux, the 
[Nill] 6.63 fim and the [Coll] 10.52 //m line fluxes. 



days it decreased considerably, and was no longer detectable af- 
ter +66 1 days. This behavior is fully consistent with the expected 
radioactive decay of 56 Co. In Figj4] we plot the calculated Co- 
masses (filled symbols) together with the predicted masses from 
the decay of 0.02 M radioactive 56 Ni synthesized in the explo- 
sion (Paper I). Since the observed Co-masses were derived from 
a single Co II transition only, these values are actually lower lim- 
its of the total Co-mass. Because 57 Co (produced by the decay 
of 57 Ni) may also contribute to the observed feature, its pres- 
ence was also included in the radioactive model, assuming solar 
and twice solar abundance ratio (dashed and dotted lines, respec- 
tively; Roche et al. 1993). We conclude that the Co-masses de- 
rived from the single [Coll] 10.52// transition agree very well 
with the model assuming that 0.02 M Q radioactive 56 Ni were 
synthesized in the explosion. 
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Fig. 6. IRAC (3.6 fim - open squares, 4.5 fim - filled squares, 
5.8 fim - open circles, 8.0 fim - filled circles), IRS PUI (open 
triangles) and MIPS 24.0 fim (filled triangles) light curves of 
SN 2004dj. The excess emission appearing on 3.6, 5.8 and 8.0 
fim peaks later at longer wavelengths, which is consistent with a 
warm, cooling dust formed one year after the core collapse. 



3.4. The continuum emission at 10 fim 

In Fig |5] the evolution of the continuum flux at 10 fim is plotted 
with the [Coll] 10.52// and [Nill] 6.63// line fluxes. There is a 
striking similarity between these curves and those of SN 1987 A 
(cf. Fig. 2 in Roche et al. 1993). The evolution of the [Nill] and 
[Co II] line fluxes were discussed above. The 10 fim flux values 
show a significant decline between +115 and 291 days, but after 
that they brighten up and by +5 10 days they reach the same level 
as in the earliest spectrum and only slightly decrease later. 

Roche et al. (1993) explained the 10 fim light curve of 
SN 1987 A as caused by optically thin free-free emission by 
~300 days, then after ~500 days the increasing flux was inter- 
preted as the sign of dust formation. We invoke the same mech- 
anism for SN 2004dj to explain the 10 fim light curve here. 
Similar temporal evolution could be identified in the IRAC and 
MIPS light curves as well, which are discussed below. 

3.5. Mid-IR light curves 

Similar to the temporal evolution of the 10 fim continuum flux 
(FigEJl, the light curves of IRAC data from the 3.6 fim, 5.8 fim 
and 8.0 fim channels also show a significant peak after +400 
days post-explosion (Fig|6]l. Such late-time MIR excess is usu- 
ally considered as a strong evidence for the presence of dust. The 
excess emission peaks later at longer wavelengths, which is con- 
sistent with a model of warm, cooling dust grains formed in the 
ejecta. Unfortunately, there are no MIPS data during the peak of 
the IRAC light curves, but the 24 fim fluxes also show a slight 
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Fig. 7. Evolution of MIR SEDs of SN 2004dj 



excess after +800 days with respect to those obtained between 
+ 100- 300 days. 

The 4.5 fim fluxes do not show the peak that the other IRAC 
channels do. The most plausible explanation for this is the pres- 
ence of the 1-0 vibrational band of CO at 4.65 pm (Section 
I3.21 i that contributes significantly to the measured flux in the 4.5 
pm channel. After — 1-500 days when the CO band disappeared 
(Fig(2j», the light curve in this channel behaves similarly to those 
in the other IRAC bands. 

In Fig|7]the evolution of the MIR SED of SN 2004dj is plot- 
ted (a vertical shift has been applied between the data for better 
visibility). This figure further illustrates the disappearance of the 
CO band as well as the shifting of the peak of the SED toward 
longer wavelengths. The fitting of SEDs with dust models is pre- 
sented in Section|4] 

We conclude that the MIR excess flux present in all IRAC 
bands and also at longer wavelengths between ~ +450 - 900 days 
is most likely caused by thermal emission of dust particles inside 
the SN ejecta. The significant increase and subsequent decrease 
of the excess emission on a timescale of a few hundred days 
argues against an IR-echo of the SN radiation from pre-existing 
CSM or ISM, because in that case the temporal evolution of the 
light curve should be at least one order of magnitude slower. 

3.6. Polarization data 

The measured true degree-of-polarization p true = 0.81 ± 0.4 % 
derived in Section 12.31 suggests a 2-sigma (weak) detection of 
polarized light from SN 2004dj in the optical (close to Z?-band) 
at +425 days. This phase is not covered well by the Spitzer ob- 
servations, but it is at the beginning of the "MIR bump"-phase, 
when the suspected dust formation just started. The detected po- 
larization, if real, fits nicely into this picture, as the photons scat- 
tered off the dust grains are expected to be polarized. 

The source of the detected polarization may be interstel- 
lar, not related to the SN ejecta. To prove or reject this hy- 
pothesis, we used the spectropolarimetric results by Leonard et 
al. (2006). They detected changing net continuum polarization 
of SN 2004dj during early phases, starting close to zero, then 
climbing up to p — 0.55 % at +90 days then declining as ~ 1/t 2 
to 0.052 % by +270 days after explosion (see FigJSJl. On the 
other hand, the measured polarization in the optical spectral lines 
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Fig. 8. Evolution of the detected true degree-of-polarization in 
SN 2004dj. Open symbols: spectropolarimetry by Leonard et. 
al. (2006); filled symbol: imaging polarimetry with HST/ ACS 
(this paper). 



was much less, going below 0. 1 % in strong lines (Ha or the Call 
triplet) and in the metallic line blends shortward of 5500 A. 

Leonard et al. (2006) also measured the interstellar polariza- 
tion in the direction of SN 2004dj, and found pism = 0.29 %. 
Obviously our detection from HST at +425 days is significantly 
higher than pjs m found by Leonard et al. Subtracting pjs m ~ 0.3 
% from the total measured polarization, psn ~ 0.5 % remains, 
which is similar to the amount of net continuum polarization de- 
tected by Leonard et al. (2006) during the photospheric phase. 

It must be emphasized that the source of the polarization dur- 
ing the photospheric phase (measured by Leonard et al. 2006) 
and the nebular phase presented here is very different. During 
the photospheric phase, the polarization is caused by Thompson- 
scattering, and the detected net polarization suggests an asym- 
metric shape of the inner ejecta (Leonard et al. 2006; Wang & 
Wheeler 2008). Leonard et al. also pointed out that this kind 
of polarization is expected to quickly decrease with time when 
the ejecta dilutes and the scattering optical depth declines, and 
this is exactly what they have observed (cf. Figj8]). On the other 
hand, the polarization detected with HS T at +425 days, when 
the ejecta are thought to be almost fully transparent and the 
electron scattering optical depth is much smaller, should be pro- 
duced by scattering on other particles. Newly formed dust seems 
to be a plausible explanation, although some kind of deviation 
from spherical symmetry is still necessary to produce observ- 
able net polarization. One possibility may be the formation of 
dust clumps distributed asymmetrically in the SN ejecta. This 
was also proposed by e.g. Tran et al. (1997) to explain the ob- 
served polarization in SN 1993J. We conclude that the detected 
polarization, psn ~ 0.5 % at +425 days is likely caused by scat- 
tering on freshly formed dust particles, which agrees with other 
signs of dust formation after this epoch. 

4. Models for dust 

In the previous section we presented several pieces of evidence 
for the presence of dust around SN 2004dj. In order to analyze 
the physical properties and estimate the total amount of dust, 
we fit theoretical SEDs from analytic and numerical models to 
the measured IRAC and MIPS data (days 267-275, 849-883, 
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1006-1016 and 1236-1246, respectively). First, we use the an- 
alytic model described by Meikle et al. (2007), then we apply 
the numerical radiative-transfer code MOCASSIN (Ercolano et 
al. 2003, 2005). Prior to fitting, the observed SED fluxes were 
dereddened using the galactic reddening law parametrized by 
Fitzpatrick & Massa (2007) assuming Ry = 3.1 and adopting 
E(B - V) = 0.1 from Paper I. The IRAC data were also cor- 
rected for the fluxes of the surrounding cluster as explained in 
Section [3~T1 We assumed a general 10 % uncertainty for all ob- 
served fluxes to represent the overall random plus systematic er- 
rors (see e.g. Kotak et al. 2005). 

4.1. Analytic models 

We fitted simple analytic models to the observed MIR SEDs 
using Eq.l in Meikle et al. (2007) assuming a homogeneous 
(constant-density) dust distribution. To estimate the dust optical 
depth (t,)), we adopted the power-law grain size distribution of 
Mathis, Rumpl & Nordsieck (1977, hereafter MRN) assuming 
m = 3.5 for the power-law index and a mm = 0.005 fim and a max 
= 0.05 fim for the minimum and maximum grain sizes, respec- 
tively. 

As there is no detectable 9.7 micron silicate feature in any of 
the observed SEDs, our initial fits using C-Si-PAH dust composi- 
tion (Weingartner & Draine 2001) failed. The lack of detectable 
SiO emission in the IRS spectra (Fig|2} also supports the absence 
of silicates in the dust composition. 

Thus, we chose amorphous carbon (AC) grains for subse- 
quent models. Values for the dust opacity ka were taken from 
Colangeli et al. (1995). For the grain density we adopted p gr = 
1 .85 g cirT 3 as the value for the AC1 material given by Rouleau 
& Martin (1991) (see e.g. Kotak et al. 2009; Botticella et al. 
2009). The grain temperature T and the grain number density 
scaling factor (k) were free parameters during the fitting. The 
dust is assumed to be distributed uniformly within a sphere, the 
radius R of which was chosen using the maximum velocity of 
the SN ejecta during the nebular phase (y max ~ 3250 km s , 
Paper II; see the method in Meikle et al. 2007), and assuming 
homologous expansion, i.e. R = v max ■ t where t is the phase of 
the SN at the moment of the SED observation. 

For SN 1987 A Wooden et al. (1993) showed that a hot com- 
ponent (T ~ 5000 K, probably caused by an optically thick gas 
within the ejecta) affects the dust-emission continuum. On the 
other hand, Meikle et al. (2007) and Kotak et al. (2009) found 
that the contribution of the hot component is very small at wave- 
lengths longer than 3 fim for SNe 2003gd and 2004et (also Type 
II-P events). In a T = 6000 K, v hot = 60 km s" 1 blackbody at 
+267 days (Kotak et al. 2009) this contribution is less than 1% 
. Late-time optical photometry of SN 2004dj was dominated by 
the flux from the cluster Sandage-96, which prevented the in- 
clusion of optical measurements in the fitted SEDs. For these 
reasons, we neglected the presence of the hot component in the 
models. 

The results are plotted in Fig. [9] Evidently a single compo- 
nent, whether blackbody or optically thin, cannot give an ad- 
equate fit simultaneously to the IRAC and MIPS fluxes. All 
single-temperature models underestimate the flux at 24 fim. 
Thus, we added a cold component to our models, similar to 
Kotak et al. (2009), which resulted in a reasonably good fit to 
all observed data. However, a single flux value obtained from 
broadband photometry may be misleading, because it may not 
be entirely caused by pure continuum emission by a blackbody. 
Since none of the IRS spectra extend to this wavelength regime, 



we could not rule out that the MIPS flux is contaminated by 
line emission. On the other hand, Kotak et al. (2009) presented 
nebular-phase MIR spectra of the Type II-P SN 2004et extend- 
ing to 30 fim that showed no significant emission lines around 
24 fim. Assuming that SNe 2004dj and 2004et had similar MIR 
spectra we considered the MIPS fluxes as caused by pure con- 
tinuum emission. 

The parameters of the best-fitting models are collected in 
Table [6] The radius of the dust-forming zone containing the 
newly formed warm dust, R warm , increased from 0.75 to 3.88 x 
10 16 cm between 270 and 1240 days, while during the same time 
its temperature decreased from 710 to 424 K. The flux from the 
cold component can be described by a simple blackbody with 
T aM = 186 - 103 K and R coM = 1.5 - 6.2 x 10 16 cm. 

The cumulative mass of freshly formed dust was calculated 
as M d = 4nR 2 T v /3K v (Lucy et al. 1989, Meikle et al. 2007) for 
each fitted epoch. The resulting masses are between 3.1 x 10~ 6 
and 1.4 x 10~ 5 M o , but it should be noted that these masses from 
analytic models are actually lower limits (Meikle et al. 2007), 
because of the implicit assumption that the dust cloud has the 
lowest possible optical depth. Theoretical studies by Kozasa et 
al. (2009) also suggested that assuming optically thin dust opac- 
ity in the MIR leads to an underestimation of the total dust mass. 
More realistic dust masses can be obtained only from numerical 
models for the optically thick dust cloud (see below). 

The last three columns in Table [6] contain the blackbody 
MIR luminosities, L warm and L co u, calculated from the radii and 
temperatures of the warm and cold components, respectively. 
For comparison, the luminosity from the radioactive 56 Co-decay 
(Leo) is also given, assuming 0.02 M Q for the initial Ni-mass 
(Paper I). Apparently the warm component dominates the MIR 
luminosity, which is comparable to Lq at ~ 270 days, but be- 
comes orders of magnitude higher than Lq by ~ 850 days. 

4.2. Numerical models 

To model the optically-thick dust cloud around SN 2004dj, 
we applied the three-dimensional Monte Carlo radiative-transfer 
code MOCASSIN. This code was originally developed for mod- 
eling the physical conditions in photoionized regions (Ercolano 
et al. 2003, 2005). The code uses a ray-tracing technique, fol- 
lowing the paths of photons emitted from a given source through 
a spherical shell containing a specified medium. For numerical 
calculations, the shell is mapped onto a Cartesian grid allowing 
light-matter interactions (absorption, re-emission and scattering 
events) and to track energy packets until they leave the shell and 
contribute to the observed SED. In newer versions it is possible 
to model environments that contain not only gas but also dust, or 
to assume pure dust regions (Ercolano et al. 2005, 2007). This 
can be applied to reconstruct dust-enriched environments of CC 
SNe and to determine physical parameters of the circumstellar 
medium (grain-size distribution, composition and geometry), as 
was shown by Sugerman et al. (2006) and Ercolano et al. (2007). 

We used version 2.02.55 of MOCASSIN for the present 
study. We chose amorphous carbon for grain composition (op- 
tical constants were taken from Hanner, 1988) for the same rea- 
son as during the fitting of analytic models described above. The 
original grid in the code allowed us to model uniform spatial dis- 
tribution of grains. Thus, we first created constant-density dust 
shells. 

We adopted four different grain-size distributions: MRN 
(with the same parameters as before) and three other cases of 
single-sized grains with radii of 0.005, 0.05 and 0.1 fim. With 
only small (r=0.005 fim) grains, we could not reproduce the ob- 
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Two components (opt. thin) 

• One component (opt. thick) 




Two components (opt. thin) 

One component (opt. thick) 




Days 849-883 




Fig. 9. Single-component blackbodies (optically thick case, dotted lines) and two-component analytic dust models (optically thin 
case, solid lines) compared to the MIR SEDs. At the first epoch, we eliminated the 4.5 fim point from the fitting, because of the 
excess flux from the 1-0 vibrational band of CO (see text for details). The IRS PUI fluxes are also shown as empty circles for 
illustration, but those were not included in the fitting. 

Table 6. Parameters for the best-fit analytic models to SN 2004dj SEDs. R/,b and Ttb belongs to the single-component blackbody, 
while the other parameters correspond to the two-component analytic model (see Section 4.1). 



Epoch 


Tbb 


Rbb 


T 

1 warm 


"warm 


Tcold 


Rcold 


M iust 


'-'warm 






(days) 


(K) 


(10 16 cm) 


(K) 


(10 16 cm) 


(K) 


(10 16 cm) 


(1(T 5 M ) 


(10 


8 erg s~ 


') 


267-275 


893 


0.19 


710 


0.75 


186 


1.5 


0.31 


19.0 


1.9 


237 


849-883 


634 


0.37 


530 


2.48 


120 


4.3 


1.11 


18.3 


2.7 


1.1 


1006-1016 


545 


0.39 


462 


2.85 


110 


4.6 


1.32 


11.1 


2.2 


0.3 


1236-1246 


490 


0.40 


424 


3.88 


103 


6.2 


1.39 


7.8 


3.1 


0.04 



Notes. Typical error is ±15 K for temperatures and ±5 x 10' 4 cm for radii. R lmn „ was calculated from v = 3250 km s 1 for every epoch (see text). 



served SEDs, but we found an adequate fitting for the other three 
cases. It agrees well with the calculations of Kozasa et al. (2009) 
that dust mass is dominated by grains with radii larger than 0.03 
/im in Type II-P SNe. 

The final results of these calculations are shown in Table [7] 
For the best-fitting models, the illuminating source was given 
as a blackbody with T BB = 7000 K and L t = 2.2-4.5 x 10 5 L o . 
The geometry of the dust cloud was assumed as a spherical shell 
with inner and outer radii R,„ and R out . For R out we adopted the 
radius of the warm component {R warm ) from the analytic models 
(corresponding to v ~ 3250 km s in velocity space). The inner 
radius was a free parameter, and the values that produced the 
best fit are shown in Table [7] 



For the dust masses we obtained values between 2-8 x 
10~ 4 M o , depending on the size and thickness of the dust cloud, 
which are more than one order of magnitude higher than the 
masses from analytic models. It should be noted, however, that in 
order to get satisfactory fitting of the MIPS fluxes, a cold black- 
body component was added to the model with the same parame- 
ters as described in Section l4~T1 

As a second step, we kept the MRN distribution, but gen- 
erated shells with p oc r~ n density profiles. We chose n — 7 
because this density distribution was found when modeling the 
optical spectra of SN 2004dj during the photospheric phase 
(Paper I). These steep density distributions in the SN ejecta are 
also predicted by numerical simulations (e.g. Chugai et al, 2007; 
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Table 7. Parameters for best-fit homogeneous MOCASSIN models to SN 2004dj SEDs 





Days 267-275 


Days 849-883 


Days 1006-1016 


Days 1236-1246 


L„(10 5 L o ) 


4.5 


4.5 


2.8 


2.2 


T BB (K) 


7000 


7000 


7000 


7000 


R in {lO [5 cm) 


3.5 


5.0 


5.0 


8.0 




2.1 


5.0 


5.7 


4.9 


Grain size: MRN 


n d , w (\0- 6 cm- 3 ) 


10.0 


10.0 


10.0 


2.0 


M A ,,„(10- 4 M Q ) 


0.2 


2.0 


4.8 


2.6 


Grain size: 0.05 fim 


>W(10 6 cm 3 ) 


10.0 


10.0 


10.0 


1.0 


M rf „ s ,(10- 4 M G ) 


0.2 


3.2 


7.6 


3.6 


Grain size: 0.1 fim 


n dus ,(\0- 6 cm" 3 ) 


2.0 


1.0 


1.0 


0.3 


M rf „ s( (10- 4 M o ) 


0.2 


2.6 


6.2 


4.2 



Utrobin, 2007). We computed symmetric grids for the density 
distribution with 10 grid points on each axis, but for compari- 
son we also produced models with higher resolution grids using 
15, 30 and 49 grid points. These models were applied for the 
SED observed between 849-883 days (see the different models 
on Fig. [TOt. The previous values for R out were again used. 

We show in Table |8]that the models with a power-law grain- 
density distribution led to somewhat different solutions, which 
resulted in dust masses one order of magnitude lower than in 
the case of constant-density models (but still 3-4 times higher 
than the results from analytic models). In order to test the effect 
of choosing n=7 on the derived dust masses, we also generated 
a model assuming p oc r~ 2 density profile (e.g. Ercolano et al. 
2007, Wesson et al. 2010). This density distribution is expected 
in a CSM produced by a stellar wind with constant outflow ve- 
locity. The resulting dust mass turned out to be a factor of 2 
higher than in the case of n=7 (see Table [8]), but lower than the 
results from the constant-density models (Table|7]i. 

4.3. Discussion of the dust models 

The results of analytic and numerical modeling presented above 
confirm that the MIR SEDs can be well explained by assum- 
ing newly-formed dust in the ejecta of SN 2004dj. The lack of 
the observed spectral features of SiO and tests of different mod- 
els yield amorphous carbon as the most likely grain composi- 
tion. The best-fitting models give large (a ~ 0.05-0.1 //m) grains, 
which is consistent with recent theoretical studies of dust in Type 
II-P SNe. Between +267 and +1246 days, the dust temperature 
decreased from ~700 to ~ 400 K. The lower limit for dust mass 
is 1.4 X 10 M , while the highest mass derived from our mod- 
els is 7.6 x 10 -4 M . The estimated amount of new dust around 
SN 2004dj is similar to those obtained for other Type II SNe, 
and therefore does not prove that CC SNe are significant dust 
sources in the Universe. 

Note that if the dust had been assumed to form optically- 
thick clumps Lucy et al. 1989), the resulting dust mass would 
have been somewhat higher. Because the mass hidden by clump- 
ing depends on the actual number and optical depth of individual 



0.1 










Analytic (x 2 = 0.73) 






Horn. MRN (x 2 = 0.80) 






Horn. 0.1 micron (x 2 = 1.17) 






Power-law distr., 1 grid points (x 2 


= 1 .09) 




Power-law distr., 30 grid points (x 2 


= 1.86) 




Power-law distr., 49 grid points (x 2 


= 2.69) 





5 10 20 30 



A. (|j.m) 

Fig. 10. Comparison of the best fitting analytic and MOCASSIN 
models with the observed SED for 849-883 days. See text and 
Tables [6j [7] and [8] for the explanation of the different models. 
Models assuming the MRN grain-size distribution were found 
to be the best to reproduce the observations. The IRS PUI flux 
(empty circle) was not used during the fitting. 



clumps, it is difficult to predict the dust mass unambiguously in 
such a model. Numerical simulations by Sugerman et al. (2006) 
and Ercolano et al. (2007) showed that the mass of hidden dust 
could be one order of magnitude higher than calculated with 
smooth density distribution. We did not investigate these mod- 
els in detail, but mention that even if the dust distribution were 
clumpy around SN 2004dj, and the dust mass were one order of 
magnitude higher than estimated above, it still would not reach 
0.01 M@. 

The modeling of the SEDs revealed a cold component as a 
T ~ 110 K temperature shell lying at 1.5 - 6 x 10 16 cm away 
from the central source. The location of this component is in 
the outer region (v ~ 6400 km s -1 ) of the SN ejecta, where the 
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Table 8. Parameters for MOCASSIN models with power-law density distribution 





L, 






Tbb 




M dusl 




(10 5 L G ) 


(10 15 cm) 




(K) 


(10~ 8 cm" 3 ) 


(10- 4 M o ) 


Days 267-275 


10 grid points 


4.5 


1.0 


7.5 


7000 


0.2 


0.1 


Days 849-883 


10 grid points 


4.0 


1.0 


25.0 


7000 


1.0 


0.4 


10 grid points 


6.0 


5.0 


5.0 


7000 


10.0 


1.3 


(n=2) 














15 grid points 


4.0 


1.0 


25.0 


7000 


3.0 


0.4 


30 grid points 


4.0 


2.5 


10.0 


7000 


800.0 


0.5 


49 grid points 


4.0 


2.0 


12.5 


7000 


300.0 


0.4 


Days 1006-1016 


10 grid points 


2.4 


1.0 


28.5 


7000 


1.0 


0.4 


Days 1236-1246 


10 grid points 


1.5 


1.0 


38.8 


7000 


1.0 


0.4 



Notes. n c i ust stands for the average number density of grains. All but one model were computed assuming n=7 power-law index in the density 
profile. The second row at days 849-883 shows the best-fit parameters derived assuming n=2 (see text for details). 




3000 AU 
1 



Fig. 11. Geometrical model of warm (inner, gray) and cold 
(outer, hatched) dust shells around SN 2004dj at ~ 850 days 

ejecta/CSM interaction is expected to take place. Such a cold 
component with similar temperature (T ~ 100 K) was found 
in SN 2004et by Kotak et al. (2009), and also in SN 2002hh, 
although the latter is slightly warmer (T ~ 300 K, Meikle et al. 
2006). However, the cold components in these Type II-P SNe 
were attributed to an IR-echo, i.e. a radiation from pre-existing 
dust reheated by the strong UV/optical radiation from the SN 
close to maximum light. 

Although pre-existing dust around SN 2004dj seems to 
be a plausible explanation for the early radio/X-ray detection 



(Chevalier et al. 2006), we do not favor the IR-echo hypothe- 
sis as the cause of the MIR cold component in this case. There 
are several arguments that do not support the probability of 
the IR-echo. First, both SNe 2002hh and 2004et were moder- 
ately or highly reddened (Ay > 1 mag), while the reddening 
of SN 2004dj was much less (Ay < 0.3 mag, Paper I). Thus, 
pre-existing CSM dust, if present, should have been much less 
dense around 2004dj than in 2002hh or 2004et. Second, the sur- 
rounding cluster, S96, is a young (f ~ 10 Myr), very compact 
cluster with a significant OB-star population (Maiz-Apellaniz 
et al. 2004, Paper II). Strong UV/optical radiation from nearby 
OB-stars should expel the ISM/dust from the cluster, although 
some CSM around massive stars resulting from mass loss via 
stellar wind might still be present. Third, the early UV/X-ray 
flash of the SN should create a dust-free cavity within a region 
of ~ 10 16 - 10 17 cm in size after core collapse and around max- 
imum light (Dwek, 1983, 1985). The radii of the dust models 
that fit the observed MIR SEDs are within this region, therefore 
the MIR radiation should come mostly from inside the cavity. 
Pre-existing dust within this region is unlikely. We also empha- 
size that the dust cloud around SN 2004et was found to have a 
roughly constant size (Kotak et al. 2009), while our models for 
SN 2004dj clearly indicate that the dust regions expand homo- 
loguously with the ejecta. For these reasons, we do not expect 
that either the cold or the warm component is caused by IR- 
echo, although some pre-existing CSM at higher distance cannot 
be ruled out completely. 

Another more probable possibility for the cold component 
is the condensation of dust grains in a CDS between the for- 
ward and reverse shocks (Section [TJ. Although this mechanism 
is supposed to work predominantly in Type Iln SNe because of 
the higher density CSM, there are hints of a similar process tak- 
ing place around SN 2004dj. Chevalier et al. (2006) pointed out 
based on radio and X-ray observations of Beswick et al. (2005) 
and Pooley & Lewin (2004) that there is a non-negligible amount 
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of CSM around SN 2004dj, probably produced by stellar wind 
from the progenitor, which has an estimated mass-loss rate of 
~10~ 6 M yr _1 . Chugai et al. (2007) presented similar results 
from studying the observed Ha profiles during the photospheric 
and early nebular phase. They argued that the observed high- 
velocity absorption 'notch' features could not form in the un- 
shocked ejecta, instead, they are supposed to be produced by the 
CDS modified by Rayleigh-Taylor instability. Applying Eq.10 
of Chugai et al. (2007), the expected radius of the continously 
growing CDS is 1.9 - 7.0 x 10 16 cm between +267 and +1246 
days in a self-similar model. These values agree well with the 
size of the cold component involved in our models (see Table 
O, which strengthens the assumption that CDS plays a role in 
producing the MIR SEDs. 

5. Summary 

Using public data of Spitzer- and Hubble Space Telescope, we 
presented a detailed analysis of MIR light curves and spectra 
on SN 2004dj between +98 and +1381 days after explosion. 
Following SN 2004et (Kotak et al. 2009), this is the second long- 
term study of the dust-formation processes around a Type II-P 
supernova. We found several pieces of evidence for dust forma- 
tion after explosion. These include 

- significant brightening in MIR light curves starting after 
+400 days, 

- detection of ~ 0.5 % polarization from the SN ejecta in the 
optical (HSTF435W filter) at +425 days. 

We fitted analytic as well as numerical models to the 
observed MIR SEDs of SN 2004dj by applying the model 
by Meikle et al. (2007) and the 3D radiative-transfer code 
MOCASSIN. The models confirmed the presence of dust as 
early as ~ +270 days and showed that the most intensive period 
of dust formation occured between ~ +270 and 1000 days. We 
found that the observed SEDs require the presence of a "warm" 
(T ~ 500 K) and a "cold" (r ~ 100 K) dust component. The 
"warm" component probably consists of freshly-formed amor- 
phous carbon grains inside the SN ejecta at v ~ 3200 km s _1 . 
The "cold" component is located at v ~ 6400 km s _1 , which 
is close to the region between the forward and reverse shock, 
where the cool dense shell is expected to be formed during the 
ejecta-CSM interaction (Chugai et al. 2007). Using smooth dust 
density distributions, the minimum dust mass was found to be 
between ~ 10~ 5 - 10 -4 M G depending on the applied model, 
but this could be one order of magnitude higher, ~ 10~ 3 M Q , 
if the dust distribution were clumpy. However, this is still sev- 
eral orders of magnitude less than what is needed to explain the 
significant dust content in star-forming galaxies. 
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